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ABSTRACT 

We report the discovery of the closest collisional ring galaxy to the Milky Way. Such rare 
systems occur due to “bulls-eye” encounters between two reasonably matched galaxies. The 
recessional velocity of about 840 km s“ ^ is low enough that it was detected in the AAO/UKST 
Survey for Galactic Ha emission. The distance is only lO.OMpc and the main galaxy shows a 
full ring of star forming knots, 6.1 kpc in diameter surrounding a quiescent disk. The smaller 
assumed “bullet” galaxy also shows vigorous star formation. The spectacular nature of the 
object had been overlooked because of its location in the Galactic plane and proximity to a 
bright star and even though it is the 60*^ brightest galaxy in the HI Parkes All Sky Survey 
(HIPASS) HI survey. 

The overall system has a physical size of ^15 kpc, a total mass of M* = 6.6 x 10^ M0 
(stars + Hi), a metallicity of [0/H]~ —0.4, and a star formation rate of 0.2-0.5 Mq yr“^, 
making it a Magellanic-type system. Collisional ring galaxies therefore extend to much lower 
galaxy masses than commonly assumed. We derive a space density for such systems of 
7 X 10“^ Mpc“^, an order of magnitude higher than previously estimated. This suggests 
Kathryn’s Wheel is the nearest such system. We present discovery images, CTIO 4-m tele¬ 
scope narrow-band follow-up images and spectroscopy for selected emission components. 
Given its proximity and modest extinction along the line of sight, this spectacular system pro¬ 
vides an ideal target for future high spatial resolution studies of such systems and for direct 
detection of its stellar populations. 

Key words: galaxies: groups: general - galaxies: interactions - galaxies: kinematics and 
dynamics - galaxies: star formation - galaxies: individual (ESO 179-13) 


1 INTRODUCTION 

Among the structures shown by interacting galaxies (Barnes & 
Hemquist 1992), the so-called collisional ring galaxies (Theys & 
Spiegel 1976, 1977) are a rare and most distinctive sub-class (Ap- 
pelton & Struck-Marcell 1996). Collisional ring galaxies occur 
from an almost head-on collision between a smaller galaxy trav¬ 
elling along the minor axis of a larger galaxy passing through, or 
close to, its centre (Appleton & Struck-Marcell 1996). The shock 
wave sweeps up and expels gas from the system, usually leaving a 
gas-poor galaxy behind (Freeman & De Vaucoleurs 1974). 

The pre-eminent example is the Cartwheel galaxy (ESO 350- 
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40; AM 0035-335) at a distance of 124 Mpc. The Cartwheel was 
discovered by Zwicky (1941), and studied in-depth by Fosbury & 
Hawarden (1977) and Amram et al. (1998). The Cartwheel shows 
a spectacular ring of star formation, with interior spokes surround¬ 
ing the host galaxy. Other well-known collisional systems include 
Arp 148 (Mayall’s Object) at 148 Mpc (Smith 1941; Burbidge, Bur- 
bidge & Prendrgast 1964; Arp 1966), IC 298 (Arp 147; Gerber, 
Lamb & Balsara 1992) at 135 Mpc, and AM 0644-741 (ESO 34- 
11 or the “Lindsay-Shapley Ring”) at 87 Mpc (Lindsay & Shapley 
1960; Graham 1974; Arp & Madore 1987). This latter object is of¬ 
ten considered to be a typical exemplar of evolved ringed galaxies 
(e.g. Higdon, Higdon & Rand 2011). 

While star-forming events resulting from galactic collisions 
may be more common than usually thought (e.g. Block et al. 2006), 
spectacular cartwheel-type galaxies, showing complete rings (but 
not necessarily with evident spokes), appear to be very rare in 
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volume-limited surveys. The previous nearest example was the Vela 
Ring galaxy (AM 1006-380 = ESO 316-43), at 72Mpc (Dennefeld, 
Lausten & Mateme 1979). Partial ring galaxies are more numerous. 
Probably the nearest in the Madore, Nelson & Petrillo (2009) cata¬ 
logue is AM 0322-374 at 20 Mpc. UGC 9893, at 11 Mpc (Kennicutt 
et al. 2008), is also listed as a collisional system, but it lacks a clear 
ring structure (Gil de Paz & Madore 2003), appearing indistinct in 
available imagery. Collisional ring galaxy systems can be identi¬ 
fied in high- resolution, deep Ha imagery where the H II regions 
of the ring are readily visible. However, Ha galaxy surveys such as 
Gil de Paz & Madore (2003), James et al. (2004), and Kennicutt et 
al. (2008) have not added to the modest number of cartwheel-type 
galaxies currently known and there are <20 convincing systems, 
e.g. Marston & Appleton 1995. Bosch et al. (2015) present the most 
recent detailed work on such systems for an apparently empty ring 
galaxy system ESO 474 G040 at ~ 78 Mpc, thought to have arisen 
from a recent merger of two disk galaxies. 

Collisional rings systems offer an excellent environment to 
study the evolution of density-wave induced starbursts (Marston & 
Appleton 1995; Romano, Mayya & Vorobyov 2008), as well as the 
dynamic interactions between the individual galaxies (e.g. Smith et 
al. 2012). A closer example of a collisional ring system would al¬ 
low its star formation, stellar distribution and evolution to be stud¬ 
ied in detail, spatially, kinematically and chromatically. Here we 
present the discovery of a new collisional ring galaxy at a distance 
of only 10 Mpc. Our paper is organised as follows. In section 2 we 
present some background and describe the discovery of the ring 
around ESO 179-13. In section 3 we present our follow-up high- 
resolution imaging and other archival multi- wavelength images. 
Optical spectroscopy is presented in section 4. Parameters for the 
system, such as masses and star formation rates, are provided in 
section 5. Section 6 presents our discussion and give our summary 
in Section 7. 


2 BACKGROUND AND DISCOVERY 
2.1 Background 

ESO 179-13 is noted as an interacting double galaxy system in 
SIMBAD, appearing to consist of a near edge-on late-type spi¬ 
ral with another more chaotic system ~ 80 arcsec to the north¬ 
east. It was first identified by Sersic (1974) and independently re¬ 
discovered by Lauberts (1982), Longmore et al. (1982) and Cor¬ 
win, de Vaucouleurs & de Vaucouleurs (1985). A detailed H-band 
image was published by Laustsen, Madsen & West (1987). 

Woudt & Kraan-Korteweg (2001) identify the most prominent 
galaxy as WKK 7460 (component A hereafter), and WKK 7463 for 
the secondary (component B). A small, third galaxy is also identi¬ 
fied as WKK 7457 to the west of component A (component C). It 
agrees with our own positional determination to 1 arcsec. The sys¬ 
tem, centred on component A, is entry 391 in the Ha catalogue of 
Kennicutt et al. (2008), but they do not list integrated Ha (- 1 - [NII]) 
flux or luminosity estimates. Lee et al. (2011), in a related paper, 
included this system in their Galaxy Evolution Explorer (GALEX) 
ultraviolet imaging survey of local volume galaxies with a quoted 
Bt = 15, but no near-UV or far-UV fluxes are quoted. Owing to 
its proximity to the Galactic plane, the galaxy was not observed by 
GALEX prior to mission completion (Bianchi 2014). 

Component A has been classified as an SB (s) m spiral in the 
RC3 catalogue (de Vaucouleurs et al. 1991) and similarly classified 
as SB (s)dm system (T = 7.5) by Buta (1995). These independent 


classifications treat the system as a distinct, late-type spiral galaxy, 
albeit in an interaction. It can be best described as a Magellanic 
dwarf or a Magellanic spiral, similar to the EMC. 

A range of heliocentric radial velocities for component A are 
available in the literature from both optical spectroscopic and Hi 
measurements. The most recent optical results are Strauss et al. 
(1992) who quote 775±36 km s“\ Di Nella et al. (1997), 750±70 
kms“^, and Woudt et al. (1999), 757±50 km s“^. These are in 
good agreement so we take an unweighted average of 761 km s“^ 
{a= 11 km s“ ^) for component A. 

The earliest Hi heliocentric velocities towards the system 
are 836±5kms“^ by Longmore et al. (1982), 840±20kms“^ by 
Tully (1988) and 843±9kms“^ by de Vaucouleurs et al. (1991). 
The Parkes Hi ‘HIPASS’ Multibeam receiver programme gives 
843±5kms (HIPASS catalogue; Meyer et al. 2004) where it 
is noted as the 60'^*' brightest in terms of integrated flux density 
(Koribalski et al. 2004). A deeper, pointed Parkes observation was 
obtained by Schroder et al. (2009). They give parameters for com¬ 
ponent A, but note that a contribution from component B is likely. 
Their heliocentric Hi radial velocity is 842 km s“ ^ from the Hi pro¬ 
file mid-point, with a velocity width of 222 km s“^ at 20% of peak. 
The integrated Hi flux density is 113.8±3.7 Jykms“^, in good 
agreement with the HIPASS value of 100. All the Hi velocities 
are consistent within small errors. These values include the com¬ 
bined contribution from all the gas in the full system. We adopt 
842±5 km s“^ for the combined Hi system velocity. 

Kennicutt et al. (2008) give a distance of 9 Mpc based on a 
Virgocentric flow model distance. We adopt a slightly further dis¬ 
tance of 10.0 ±1.0 Mpc, based on Kennicutt et al. (2008), but scal¬ 
ing from a Hubble constant Hq of 75 to 68 km s“^ Mpc“^ due to 
the latest WMAP (Hinshaw et al. 2013) and Planck results (Ade 
et al. 2014). At this distance, 1 arcsec equates to 48.5 pc. Wouldt 
& Kraan-Korteweg (2001) give a modest colour excess of i?(B- 
V)= 0.25 mag in the surrounding field, giving Ay = 0.78 mag 
while Schlafly & Einkbeiner (2011) give Ay ~ 0.75 mag for 
Galactic extinction in this general direction, which we adopt here¬ 
after. We will also use Ahq = 0.61 mag based on the reddening 
law from Howarth (1983). 

Apart from the basic information above, the system has been 
little studied because of its location in a crowded, low-latitude 
area near the Galactic plane (|6| =-8 degrees) and proximity to 
the bright, 7.7 mag AO IV star, HD 150915, located 72 arcsec due 
south. There are actually two bright stars here, but only HD 150915 
is listed in SIMBAD. The second star is ~ 10 arcsec north of 
HD 150915 and ~2 magnitudes fainter. It was undoubtedly blended 
with HD 150915 in the original photographic imagery. 

2.2 Discovery of the Ring 

The low heliocentric velocity meant the system was discovered as 
a strong Ha emitter in the Anglo-Australian Observatory (AAO) / 
U.K. Schmidt Telescope (UKST) SuperCOSMOS Ha survey of the 
Southern Galactic Plane (SHS; Parker et al. 2005) during searches 
for Galactic planetary nebulae (Parker et al. 2006). The bandpass 
of the survey filter (Parker & Bland-Hawthorn 1998) includes all 
velocities of Galactic gas and emission from extragalactic objects 
out to the Virgo and Fornax galaxy clusters at ~1200kms“^. 

Fig.[2shows the 4x4 arcmin discovery images from the SHS 
survey data centred on the system. From left to right there is the 
Ha image, the matching broadband short-red ‘SR’ image and the 
image obtained by taking the SR image from the Ha image (ef¬ 
fectively a continuum subtracted image). Component A is at the 
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Figure 1. Left panel: The 4x4 arcminute SHS broad-band short-red ‘SR’ image centred on component A at 16h47ml9.7s, -57°26’3r’ (J2000). North is top 
and east left. Mid panel: The matching Ha image. Right panel: The continuum subtracted image obtained by taking the SR from the SHS Ha (-H [Nil]) image 
scaled at 99.5%. North is at the top and east to the left. Component B is labelled to the north-east and component C is indicated with an arrow. 


centre. Component B is labelled to the north-east while the smaller 
compact component C is indicated with an arrow. 

The SHS continuum subtracted image reveals a spectacular 
ring of emission knots not obvious, without context, even in the 
Ha image while the main component A galaxy has almost disap¬ 
peared. The ring is approximately (but not perfectly) centred on 
component A. Surprisingly, given its previous spiral morphologi¬ 
cal classification, component A does not show up well in the SHS 
continuum subtracted image and appears to show little star forma¬ 
tion. This indicates that along with its apparently highly elongated 
nature (so not an elliptical) it is has been largely stripped of gas. 
The smaller, more irregular galaxy (component B) shows exten¬ 
sive, clumpy Ha emission. A third, faint compact system (compo¬ 
nent C) is seen 38 arcsec to the west of component A, identified as 
WKK 7457. This object is clear in the broad-band B, R and /-band 
SuperCOSMOS images and is also undergoing star-formation. 

The strong Hi detection noted earlier indicates that a signifi¬ 
cant, neutral Hi component still remains within the overall system 
environment. The lack of much star formation in component A sug¬ 
gests that most of this Hi may be located outside of its inner disk. 


3 FOLLOW-UP IMAGING AND MULTI-WAVELENGTH 
COMPARISONS 

Discovery of this collisional ring system led us first to investigate 
and compile existing multi-wavelength data from the archives. We 
then obtained deeper, higher resolution imaging and spectroscopy 
of the major emission components (see below), to provide estimates 
of some key physical parameters for this system. 

3.1 CTIO 4iti-telescope Imaging 

High resolution imagery at the Blanco 4-m telescope at CTIO was 
obtained in June 2008 using the wide field mosaic camera. We used 
five filters: U-band, Ha, Ha off-band (80A redward of Ha), [O III] 
, and a (wider) [Olll] off-band filter centred at 5300A. The field 
of view is 30 arcmin on a side and the plate scale 0.26 arcsec px“ ^. 
The seeing measured from the data frames was typically 1.2 arcsec. 
Exposures were 5 minutes for the emission-line and off-band filters, 
and 60s for the U-band filter. The airmass during the exposures 


averaged 1.3. To subtract the continuum, the Ha off-band image 
was scaled down by a factor of 0.88, and the [O III] off-band by 
a factor of 0.14, to compensate for the filter widths. These factors 
were established by minimizing residuals for a set of selected field 
stars. For [O III] that works very well and results are in good agree¬ 
ment with the filter curve. For Ha this can be an issue if the chosen 
stars have Ha in absorption. This does not appear to be the case 
for the field stars chosen which provide a consistent scale-factor 
to ~ 2%. It is best to minimize stellar residuals rather than filter 
curves as stars are a more important contribution to the continuum 
than bound-free emission from Hi regions. 

A montage of these CTIO images is shown in Fig. They 
reveal the distribution and variety of the emission structures as well 
as data for quantitative flux estimates. From these deeper, higher 
resolution CTIO Ha and [O III] continuum subtracted images the 
distribution of ionised gas is seen to be far more complex than is 
evident from the SHS discovery data. 

The leftmost panel is a 3 x 3 arcminute CTIO Ha image with 
the off-band frame subtracted. We avoid the bright star to the south 
where CCD blooming is serious. Emission associated with com¬ 
ponent C and at the southern tip of component A’s disk previously 
gave the impression of a circular ring in the SHS as the bright 
star obscured emission further south. The higher resolution CTIO 
data reveals an elliptical emission ring with a major-axis diameter 
of 127 arcsec (~6.2kpc). The inscribed oval was positioned to fit 
the prominent emission now seen to extend further south. The el¬ 
lipse is not centred on component A, but at 16h47ml9.5s, —57°26’ 
44” (J2000), ~ 13 arcsec (or ~0.6 kpc) south-south-east. Some low- 
level Ha emission is now seen across component A, surrounding 
component B and with faint, localised emission features and blobs 
around the entire system. An interesting feature is the elongated na¬ 
ture of some of these emission knots in a north-east direction from 
A to B. The main [Olll] emission features follow the equivalent 
Ha structures though at typically a third to a half of their native 
integrated pixel intensities. 

The mid panel is the matching CTIO continuum subtracted 
[O III] image (narrow on and off-band [O III] images used). Both 
these images are presented at 90% linear scaling to reveal the full 
extent of the emission (this sets the upper and lower limits based 
on the 90% pixel intensity level where a histogram of the data is 
created and the limits are set to display the percentage about the 
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mean value). With full pixel range, saturated pixels from the 7.7th 
magnitude star HD 150915 would prevent detail being seen. Low 
level diffuse [O III] flux is seen for component A but it is not con¬ 
spicuous. This shows it has almost been completely cancelled out 
indicating component A is mostly composed of normal starlight. 

The right hand panel is the off-band [Olll] filter image 
(5300A) that effectively represents the H-band starlight without 
the [O III] -rH/3 emission line contribution. This shows dust lanes 
along and at the south-west extremity of the disk, concentrated 
starlight from components A, B and C and an envelope of diffuse 
starlight around the entire system. This “common envelope” was 
already noted by Laustsen et al. (1987) from old 73-band photo¬ 
graphic data taken with the ESO 3.6-m telescope in the 1980s. This 
raises an interesting question and also a possible explanation for 
why the apparent axial ratios of the highly elongated component A 
differs so much from the surrounding oval ring. Currently known 
collisional ring systems with low impact parameter have axial ratios 
that are broadly similar to that of the target galaxy. Such arrange¬ 
ments are a natural consequence of the copious gas and ISM being 
swept up and compressed by the resultant density wave having a 
similar radial (and perhaps also vertical) distribution as the target 
galaxy stellar content. This does not appear to be the case here un¬ 
less what we are actually seeing is the residual, less-inclined bar of 
a more face-on disk system that extends out to the entire ring and 
whose presence can still be seen in the low-level diffuse starlight 
refereed to earlier. If component A is not a residual bar but a distinct 
late-type disk galaxy (as commonly assumed) stripped of gas, then 
the gas distribution could have become less flattened somehow. Al¬ 
ternatively the original target galaxy could have been stretched out 
along the direction of the impact towards component B before, dur¬ 
ing or even after the bulls-eye collision occurred. Though compo¬ 
nent B is only 13% of the mass (in stars) of component A these 
effects could have been strengthened due to the low ‘Magellanic’- 
type masses of the system as a whole. 

Finally, Fig. is a composite Red, Green Blue (RGB) CTIO 
image displaying the system’s spectacular nature. The red channel 
is Ha, the green [O III] and the blue broadband V. Component A 
is mostly star-light with dust lanes evident in this higher resolu¬ 
tion data. The few small Hll knots appear red, likely due to dust 
extinction within the disk. 

Based on the SHS data, literature data and our new CTIO im¬ 
ages, Table[^presents some key parameters for galaxy components 
A, B and C, including the central position of the best fit ellipse to 
the main emission ring. The position for component A is tied to 
the 2MASS (Skrutskie et al. 2006) near-IR imagery and astrome¬ 
try. The positions refer to the centre of each assumed galaxy. The 
centre of component A is derived from the dominant starlight in the 
7?-band. The separation between the centres of components A and 
B is ~80arcsec. Using the distance of lOMpc gives a projected 
separation between them of 3.9 kpc. 

Component A has a major-axis of ~86arcsec (4.2 kpc) mea¬ 
sured to the outer edge in the stellar continuum from off-band 
[Olll] and U-band CTIO images where a decline in the back¬ 
ground is seen before the star-forming ring is encountered. When 
checked against the broad band R, I and 2MASS images when au¬ 
tocut at the 90% pixel intensity level the disk appears to extend 
only to 70-80 arcsec. At low isophotal thresholds starlight extends 
out to the star-forming ring both along the major and minor axes. 
Component B has a major axis of ~33 arcsec (1.6 kpc) determined 
in the same way. The strongest Ha emission in B is also in an oval 
structure of nucleated knots ~32 arcsec across. This is in excellent 
agreement with the value estimated from the concentrated starlight. 


Starlight internal to this emission oval is clear in Fig.j^and B may 
extend to 70 arcsec (~3.4kpc) in low level, isolated clumps of star 
formation and diffuse starlight. 

The overall system envelope, including all assumed associ¬ 
ated broad and narrow band emission and diffuse starlight as seen 
across the region in the right panel of Fig.|^ has a major diameter 
of 318 arcsec and a minor axis of 198 arcsec, giving a total system 
size of ~15.4x9.6kpc. 


3.2 Near and mid infrared imaging 

The 2MASS colour combined J, H and Ks band image (as red, 
green, blue) is presented in the left panel of Fig. The old stel¬ 
lar population of component A is clear in the near-IR and seen to 
be well separated from component B (surrounded by a blue circle) 
located towards the top left of the image, indicating it is a dis¬ 
tinct galaxy. Component C (surrounded by a blue circle to the right 
of Component A) is also just visible. Each near-IR component is 
scaled at the 95% level in the linear combination. A low level back¬ 
ground is also apparent which needs to be included in any overall 
near-IR magnitude estimates. The NIR magnitude differences be¬ 
tween components A and B is A J ~ ATT ~3.1. There is also 
only a mild colour term for each main component with J—K~0.3 
for component A and —0.05±0.05 for component B. 

The Wide-field Infrared Survey Explorer, WISE (Wright et 
al. 2010) surveyed the entire sky at 3.4, 4.6, 12 and 22/rm (bands 
named W1 to W4). The right panel of Fig. shows the WISE ‘hi¬ 
res’ mid-IR combined colour image of bands W4, W3 and W2 (as 
blue, green, red) from drizzled data. The WISE mission provides 
“Atlas” images, co-added from individual frames, through the web 
server hosted at IPAC[^ The images have beam sizes of 8.1, 8.8, 
11.0, and 17.5 arcsec for the Wl, W2, W3 and W4 bands, respec¬ 
tively. 

To regain valuable spatial resolution we made use of the 
ICORE (Image Co-addition with Optional Resolution Enhance¬ 
ment software (Masci 2013) which exploits the stable point 
spread function (PSF) of WISE. The “drizzled” co-added images 
make use of a VPLR (Variable Pixel Linear Reconstruction) tech¬ 
nique employing a “Tophat” PRF (Point Response Function) ker¬ 
nel. The “HiRes” images additionally use a MCM (Maximum Cor¬ 
relation Method) technique as outlined in Masci & Fowler (2009). 
The drizzle images have a resolution of 5.9, 6.5, 7.0 and 12.4 arcsec 
for bands Wl to W4 respectively (Jarrett et al. 2012). Component A 
is prominent in W3 centred at llfim but with sharp cut-on and cut¬ 
offs at 7 and 17 fim. This is likely indicative of thermal emission 
from diffuse interstellar dust and polycyclic aromatic hydrocarbons 
(PAHs) in the disk. The presence of clear dust lanes visible in the 
CTIO U-band image along component A supports this and shows 
that component A has not been entirely stripped of its gaseous com¬ 
ponent. 


3.3 Existing radio data 

In the radio regime the system is detected in the second epoch Mo- 
longlo Galactic Plane Survey (MGPS2; Murphy et al. 2007). The 
MGPS2 was performed with the Molonglo Observatory Synthesis 
Telescope at a frequency of 843 MHz and with a restoring beam of 


^ http://irsa.ipac.caltech.edu/applications/wise/ 

^ http://web.ipac.caltech.edu/staff/ fmasci/home/icore.html 
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Figure 2. Selected 3x3 arcmin regions from the CTIO 4m imaging centred at 16h47m21s, -57°26’03” (J2000). Left panel: the CTIO continuum subtracted 
Ha. The green oval has a major axis diameter of 127 arcsec (~6.2kpc) and a minor axis of 91 arcsec (~4.4kpc). Mid panel: matching CTIO continuum 
subtracted [O III] image. Both these images are at 90% scaling to better reveal the emission otherwise the 7.7th magnitude star dominates. Right panel: CTIO 
[OlII] off-band image. This represents the V-band starlight without any [OlII] and H/3 emission line contribution. Dust lanes can be clearly seen at the 
southern tip of component A. The more diffusely distributed starlight can also be seen in and around components A and B. 


Table 1. Fundamental parameters of Kathryn’s wheel. Component A is the main system near the ring’s centre. Component B is the smaller system directly to 
the north-east. Component C is projected internal to the ring to the south-west. It has also has a concentrated stellar component so is also considered a distinct 
galaxy. The oval ring completely encompasses component A. 


Parameter 

A 

B 

C 

Ring 

A-B separation 

RA (J2000) 

16 47 19.7 

16 47 27.8 

16 47 15.5 

16 47 21.1 


DEC (J2000) 

-57 26 31 

-57 25 35 

-57 26 44 

-57 26 36 


Size (arcsec) 

86 

33 

16 

127 

80 

Size (kpc) 

4.2 

1.6 

0.8 

6.2 

3.9 


~45 arcsec. In MGPS2, the system is resolved into two main com¬ 
ponents catalogued as MGPS J164718-572710 (associated with 
component A) and MGPS J164728-572544 (components). The re¬ 
ported peak and integrated flux densities are 15.4±1.1 mJy/beam 
and 62±6 mjy for the component A and 11.5±1.1 mJy/beam and 
34±5 mJy for components. These radio counterparts are partially 
resolved, as seen in Fig.|^ Components shows a small north-east 
extension which doesn’t correspond to either the optical or MIR 
emission. Component A peaks in the region obscured by the bright 
foreground star HD 150915 and also shows two fainter extensions 
corresponding to the emission in both the optical and MIR (Fig. 

left and middle). No significant radio emission at 843 MHz is 
observed from the inner region of the star forming ring associated 
with component A. It is possible that component C also has a radio 
contribution hinted at by the overlaying elongated MGPS2 radio 
contours in this region. 

The system is also detected in the Parkes-MIT-NRAO (PMN) 
South Survey (PMN; Wright, Griffith, Burke and Ekers 1994) at 
4.85 GHz, and catalogued as PMN J1647-5726. Because of the 
much poorer PMN resolution (~4.2 arcmin) the system is only 
detected as an unresolved point source with a flux density of 
37±8 mJy, very close to the detection limit of 32 mJy. Fig.[^(right) 
shows the PMN radio contours with a peak flux of 37 Jy/beam. An 
apparent extension to north-east, similar to the one at 0.843 GHz, 
is visible. Due to the large beam size of the PMN survey, it is likely 
that this extension is caused by the confusion with the nearby radio 
point source PMN J1648-5724 ~5.5 arcmin to the north-east. This 
is intriguing as the MGPS2 contours at the north-east edge of com¬ 


ponent B and the south-west edge around the nominal radio point 
source PMN J1648-5724 are almost aligned. However, there is no 
evidence of an optical galaxy counterpart to this compact source. 


Using the integrated flux densities at 4.8 GHz and 0.843 GHz 
(summed components A and B) we estimate a system-wide spec¬ 
tral index of a ~-0.5 (defined from Sv oc v°‘). This estimate is 
uncertain for several reasons: (i) the accuracy of the flux density 
at 4.8 GHz is poor because of possible confusion with the nearby 
source and proximity to the sensitivity limit, (ii) the accuracy of the 
flux density at 0.843 GHz suffers from the measurement method 
(Gaussian fitting) which assumes that both radio objects are unre¬ 
solved and (iii) the estimate of the spectral index from only two 
points assumes that the spectral energy distribution (SED) is flat in 
that spectral region (e.g. the free-free emission in the v < 1 GHz 
region could be optically thick). However, the estimated radio SED 
is in agreement with SED’s of normal galaxies (Condon 1992); i.e. 
a mixture of thermal emission originating from bright, star form¬ 
ing, Hll regions, and non-thermal emission from diffused, ultra- 
relativistic electrons and supernova remnants (SNRs). 


A proper interpretation of the components’s radio SED and 
separation of the thermal and non-thermal component (which could 
lead to an independent estimate of the SFR) requires new, high res¬ 
olution and high sensitivity observations at several frequencies and 
in the high frequency range (1-10 GHz). 
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Figure 3. A colour composite RGB CTIO image displaying the system’s spectacular nature. The red channel is H«, the green [O III] and the blue is broadband 
V. Component A is mostly normal staidight with dust lanes also evident in this higher resolution data as well as the likely effect of dust on the small Hll knots 
along the disk. 


4 FOLLOW-UP OPTICAL SPECTROSCOPY 

Initial spectroscopy was obtained using the now decommissioned 
Double Beam Spectrograph (DBS) on the 2.3 m MSSSO telescope 
in June 2005. The four positions were centred on the brightest 
knots in the system and are summarised in Table and shown in 
Fig|^ The DBS gives two simultaneous blue and red spectra and 
our choice of gratings covered the range 3520-5430A in the hlue 
and 5950-6900A in the red. The dispersion was 1 A per pixel in the 
hlue arm and 0.5 A per pixel in the red arm. The 1-D spectra were 
obtained by summing over the extent of the compact emission on 
the slit and were flux-calibrated following standard IRAF routines. 
The typical on-site seeing is 1.5-2.5 arcsec. 

Additional spectroscopy was obtained using the WiFes 
integral-field instrument on the 2.3m MSSSO telescope (Dopita et 
al. 2010) in July 2011. WiFes is an ‘image slicing’ integral-field 
spectrograph made of up 25 1 arcsecwide adjacent ‘slices’ on the 
sky each of which is 38 arcsec long. This creates an effective on-sky 


field of view of of 25 x 38 arcsec with the long axis oriented east- 
west in standard configuration. The sky background was subtracted 
using empty areas of the observed field or from offset sky point¬ 
ings. The weather was partly cloudy and though flux calibration 
was performed, the atmospheric transmission changed significantly 
between pointings. We used the R = 7000 grating with coverage 
from 5800A to 7000A and dispersion of 0.44A per pixel. Only the 
prominent knots in each WiFeS pointing were extracted and used 
as the S/N in the inter-knot regions is too low (partially because of 
the non-optimal observing conditions) for additional useful spectra 
to be extracted or to create a 2-D line intensity map or a reliable 
2-D velocity field. Deeper Integral Field Unit (IFU) observations 
on a larger telescope are required. 

A total of seven combined DBS and WiFeS pointings were 
obtained, centred on the brighter knots in the system as shown in 
Fig. 0 including three WiFeS repeats of some of the earlier DBS 
observations for consistency checks. The positions are listed in Ta¬ 
ble]^ For WiFeS several knots can be observed simultaneously for 
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Figure 4. Left: a 4x4 arcmin colour mosaic of the 2MASS J, H and Ks band data of the system (as blue, green, red). North is top and east left. Component A 
at the centre is clearly visible revealing the old residual stellar population. Component B (surrounded by a blue circle) is obvious towards the top left of the 
image, indicating it is a distinct galaxy. Component C (surrounded by a blue circle to the right of component A) is also just visible. Right: The equivalent 
‘hi-res’ WISE mid- IR combined colour image of bands 4, 3 and 2 (as blue, green, red) from drizzled data. Component A is strong in W3 centred at 12fim. 
The emission is likely thermal from diffuse interstellar dust and PAHs. Several knots in oval ring are particularly bright in W4. This is related to hot (~100 K) 
dust continuum emission from ongoing star formation in these regions 



Figure 5. Left and middle: SHS Ho: quotient image and WISE RGB colour composite, respectively, overlaid with contours from the MGPS2 radio image 
at 3, 5, 7, 9, 12 and 15 X(7mgps2 {o'mgps 2 = 1-1 mJy/beam). These panels are 5x6 arcmin in size. Right: MGPS2 total intensity image overlaid with 
radio contours from MGPS2 (same as in previous images) and from PMN at 3, 4 and 5 Xapmn {(Tprun = 8 mJy/beam). Here the image is 13x18 arcmin. 
North-east is to the to left in all panels. 
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each pointing. A spectrum through the data-cube was obtained by 
placing apertures over each knot and wavelength summing the pix¬ 
els within each aperture. A combined 1-D spectrum was produced 
for each knot and used to derive velocities from fits to the [N II], 
Ha and [S II] emission lines and taking the weighted average. 

Velocities to better than ~50kms“^ are achievable with the 
DBS system and grating combination, but with some flop and flex¬ 
ure being an issue with this large rotating spectrograph off the Nas¬ 
myth focus (Rodgers et al. 1988). Spectral calibration errors from 
the red WiFeS datacube are of the order of O.OSA (see Childress 
et al. 2014) or 2.3kms“^ at Ha. Arc calibrations exposures were 
taken for each object observation to remove any effect of tempera¬ 
ture variations inside WiFeS. This set-up enables us to resolve ve¬ 
locity dispersions across resolved systems down to a ~ 20 kms“^ 
(Vogt et al. 2015). In our data velocity measurements of individual 
knots in one WiFes pointing agree to better than 10kms“^. DBS 
velocities show significant scatter of 50kms between the blue 
and the red arm and are less accurate. Where the same knot is mea¬ 
sured with both WiFes and DBS there is acceptable agreement in 
two cases but the offset is fairly large (~50km s“^ in the third case 
(Knot IV) and is at the upper error level that might be expected. It 
is clear from Fig. |7]that the spectral continuum for Knot IV is lower 
than for the other spectra shown. 

Fig. [^presents the 1-D blue and red emission line spectra from 
the DBS for knots II, III, IV and V. Extracted regions are shown 
enlarged to show the fainter lines. The spectra are typical of Hll 
regions and exhibit the standard emission lines of [O II] 3727A, the 
Balmer series (HJ H 7 , H/3) and [Olll] in the blue region. In the 
red there is Ha and the weak bracketing [N II] lines. The [S II] 
doublet is also well detected and resolved providing an opportu¬ 
nity to estimate electron densities. These are near the low density 
limit (rie ~2-100 cm“^) from measurements across four knots. For 
knot III, which is the highest surface brightness compact emission 
region in the entire system, the higher S/N spectrum reveals [Ne III] 
lines at 3869 and 3968A rest wavelengths. There is no evidence for 
Wolf-Rayet features in any of these H II region spectra. 



35.0 30.0 25.0 16:47:20.0 15.0 10.0 05.0 

RA J2000 

Figure 6. The locations of the spectroscopically observed emission knots 
given in Table|^overlaid on the 4x4arcminute SHS Ho quotient image. 


Table 2. Summary DBS and WiFeS IFU observations from June 2005 
(DBS) and July 2011 (WiFeS). The DBS velocities are accurate to ~ 
50 km s“^ and WiFeS ~ 20 km s“^ 


Pointing 

RA (J2000) 

DEC (J2000) 

Instrument 

Velocity 

Knot I 

I.l 

16:47:28.68 

-57:25:34 

WiFes 

kms“^ 

799 

1.2 

16:47:26.61 

-57:25:34 

WiFes 

812 

1.3 

16:47:27.15 

-57:25:29 

WiFes 

808 

1.4 

16:47:26.60 

-57:25:21 

WiFes 

807 

Knot II 

II. 1 

16:47:27.90 

-57:25:45 

DBS 

803 

II. 1 

16:47:27.90 

-57:25:45 

WiFes 

799 

II.2 

16:47:28.90 

-57:25:33 

WiFes 

800 

Knot III 

III.l 

16:47:25.30 

-57:26:04 

DBS 

806 

III.l 

16:47:25.30 

-57:26:04 

WiFes 

823 

III.2 

16:47:25.80 

-57:26:15 

WiFes 

836 

III.3 

16:47:27.10 

-57:26:07 

WiFes 

834 

Knot IV 

IV. I 

16:47:19.68 

-57:25:51 

DBS 

827 

IV. I 

16:47:19.68 

-57:25:51 

WiFes 

880 

IV.2 

16:47:21.50 

-57:25:48 

WiFes 

864 

IV.3 

16:47:20.90 

-57:25:49 

WiFes 

871 

KnotV 

VI 

16:47:17.18 

-57:26:07 

WiFes 

909 

V.2 

16:47:17.82 

-57:26:03 

WiFes 

901 

Knot VI 

VI. 1 

16:47:14.80 

-57:26:23 

DBS 

751 

Knot VII 

VII. 1 

16:47:15.90 

-57:26:44 

WiFeS 

938 

VII.2 

16:47:18.25 

-57:26:51 

WiFes 

910 

VII.3 

16:47:14.83 

-57:26:48 

WiFes 

959 

VII.4 

16:47:17.55 

-57:26:55 

WiFes 

932 


5 ANALYSIS 

5.1 Optical Photometry 

Component A (the main ESO 179-13 galaxy) has optical photo¬ 
graphic photometry in the Bj, R and /-bands (Doyle et al, 2005) 
while Lee et al. (2011) give B ~15. In the broad-band SuperCOS- 
MOS photographic data (Hambly et al. 2001a) the isophotal im¬ 
age parameters for component A are blended with the bright star in 
the /Jj-band so no reliable magnitude is possible from the stan¬ 
dard image-analysis mode (lAM) processing. This led Doyle et 
al. (2005) to estimate individual Sextractor magnitudes for com¬ 
ponent A when constructing optical photometry for HIPASS galax¬ 
ies. These were bootstrapped to the various SuperCOSMOS zero- 
points for each band. However, the areal profile fitting and image 
deblending techniques employed by the SuperCOSMOS lAM soft¬ 
ware (Hambly et al. 2001a) do enable component magnitude esti¬ 
mates to be obtained. There are three //-band magnitudes available 
from the first and second epoch broad-band surveys and from the 
SR exposure taken as part of the Ha survey, /-band SuperCOS¬ 
MOS magnitudes are also available. For isolated sources such mag¬ 
nitudes are accurate to ±0.1. These results are presented in Table|^ 
The lAM results from components B and C are likely to be more re¬ 
liable than for component A whose /-band magnitude in particular 
seems grossly underestimated when compared to the later 2MASS 
near- IR photometry and any sensible colour term. 
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Figure 7. The four DBS spectra, numbered according to position in Table|^ Extracted spectral regions are shown enlarged to bring out fainter lines. The red 
and blue zoom windows show the spectral ranges 3700^400A and 6500-7700A. The spectra are typical of Hll regions with strong Balmer emission and 
forbidden lines of [O III] , [Nil] and [S11]. Note the [Neill] lines in the blue data of the higher S/N spectrum for knot Ill. 


5.2 Infrared Photometry 

The available infrared photometry obtained is presented in Table 
l^with data gathered from various archives including Skrutskie et 
al. (2006) for 2MASS. The integrated Tfs-band photometry was 
obtained from the 2MASS ‘Large Galaxy Atlas’ (LGA: Jarrett et 
al. (2003). A total It's-band magnitude of 8.95 is obtained from 
an aperture of 177.8 x 56.9 arcsec that combines the three sys¬ 


tems though the contribution from component C and diffuse com¬ 
ponent is minimal. We carried out differential aperture photometry 
on the two main sub-systems A and B and found a TTs-band dif¬ 
ference of 2.2 mag with component A being the brighter. We find 
m_R'(A) = 10.3 mag and = 12.5 mag, for the main bod¬ 

ies of components A and B only (i.e. the main disk for A and the 
core region of B). 

The system corresponds to IRAS 16430—5721 in the Point 
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Table 3. SuperCOSMOS and Sextractor isophotal magnitude estimates 
for main system components. The asterisked I band entry for compo¬ 
nent A came from the lAM data from the SR/Ha SuperCOSMOS images. 
^ Sextractor values ai'e from Doyle et al. (2005). 


Component 

R1 

R2 

I 

Origin 

A 13.1 

- 

13.3 

12.8 

Sextractor^ 

A 

14.98 

13.85 

11.3* 

lAM (SuperCOS) 

B 

15.61 

15.26 

13.86* 

lAM (SuperCOS) 

C 

16.27 

16.26 

16.0 

lAM (SuperCOS) 


Source Catalog for the Infrared Astronomical Satellite (Beichman 
et al 1988). The galaxies were detected in the 60 and lOOp-m bands, 
but only upper limits are reported in the 12 and 25 pm bands. 

Both the AKARI infrared camera (IRC) All-Sky Survey Point 
Source Catalogue (Ishihara et al. 2010) and the AKARI far infrared 
surveyor (FIS) Bright Source Catalogue (Yamamura et al. 2009) 
also report flux densities for the objects. However, the AKARI 9p.m 
detection comes from the bright star and not the galaxy. HD 150915 
is a visual binary with components lOarcsec apart. The companion 
star (not in SIMBAD) is redder with (J-Ks) = 1.1. The two stars 
have similar magnitude in the mid-infrared: 7.45 and 7.36 in W2, 
and 7.20 and 7.22 in W3. The AKARI 9p.m position is in between 
the two stars. The two stars have no infrared excess (McDonald, 
Zijlstra & Boyer 2012). The AKARI FIS measurements at 65 and 
I60pm are < 3cr detections. 

The Planck Catalog of Compact Sources Release 1 (PLANCK 
collaboration: Ade et al. 2014) lists flux densities at 545 and 
857 GHz (550 and 350pm ), although the 545 GHz measurement 
is only a ~2cr detection. A summary of the available infrared pho¬ 
tometry is presented in Table|^ 

5.3 Hi and optical spectroscopic velocities 

The Hi HIPASS spectrum of the system from a pointing at 
16h47m20.0s, -57°27’06” (J2000) is reproduced in Fig. [S] The 

Table 4. Infrared photometry for the combined system unless othei'wise 
stated. WISE values are aperture conected. 


A 

Band 

Flux density [mJy] 

[mag] 

1.2|j.m 

2MASS I 


9.571 

1.7 pm 

2MASS H 


9.122 

2.2pm 

2MASS Ks LGA 

184 

8.895 

3.4pm 

WISE W1 

129.51 

8.45 

4.6|j.m 

WISE W2 

82.38 

8.30 

12|j.m 

WISE W3 

269.80 

5.17 

22|j.m 

WISE W4 

739.20 

2.63 

12|j.m 

IRAS 12|j.m 

< 250 


25 M-m 

IRAS 25 M-m 

< 295 


60 pm 

IRAS 60|J.m 

3930 ± 830 


100 pm 

IRAS 100pm 

14200 ± 1600 


65 pm 

AKARI FIS N60 

2179 ± 748 


90 pm 

AKARI FIS WIDE-S 

4982 ± 272 


140pm 

AKARI FIS WIDE-L 

8730 ± 1140 


160 pm 

AKARI FIS N160 

6105 ± 3460 


350p.m 

Planck 857 GHz 

7264 ± 1629 


550p.m 

Planck 545 GHz 

1512 ± 760 



20 % line width is 222 kms“^ centred at a systemic heliocenPic 
velocity of 842km s“^ (Koribalski et al. 2004). 

For our optical observations we used the WiFeS and DBS 
spectra to derive heliocentric velocities of several of the most 
prominent emission knots, using weighted averages from fits to the 
key emission lines. These results can be used to trace the velocity 
of the knots in the main ring and for the other star forming regions 
in component B, but not of the central disk where there is little star 
formation and where we did not obtain spectra. We rely on the pub¬ 
lished results for component A. 

Our new spectroscopic results are shown in Fig.|^ The verti¬ 
cal axis shows the heliocentric velocity of each knot. On the RHS 
of the figure we have added the HIPASS spectrum plotted with ve¬ 
locity on the same scale for ease of comparison. The horizontal 
axis shows the angle along the ring, where a position angle of zero 
corresponds to the direction of galaxy component B. This is consid¬ 
ered the bullet because of the lack of other obvious galaxies in the 
broad-band imaging within 30arcmin (~5.2Mpc). Black symbols 
are knots which we consider as part of the ring as they fall along 
the ring’s locus even if they have a discrepant velocity. Red symbols 
are for knots in galaxy components B or C, located inside or out¬ 
side of the ring. Circles are the WiFeS data and triangles the (less 
accurate) DBS data. There appears to be one particularly anoma¬ 
lous DBS velocity for knot VI at only 750km s“^. This knot is the 
second most intense star-forming region in the entire system. There 
is no corroborating WiFeS spectra for this knot to check the DBS 
result. 

The velocity distribution for the knot samples around the ring 
shows a sinusoidal pattern as function of ring position angle, with 
systemic velocity of 880 kms and an amplitude of 60kms“^. 
This velocity pattern can interpreted as rotation or expansion, or 
a combination of both. The fact that the extrema of the velocity 
fall along the major axis is indicative of rotation being the domi¬ 
nant component. The velocity corrected for inclination of the ring 
gives 85kms“^. If this is rotation, the interior mass would be 
~ 10^° M 0 , a little higher than derived above for component A. 

The red points at angle zero show the velocity distribution of 
the most conspicuous knots seen in the component B galaxy (the 



500 1000 


Figure 8. The Hi spectmm from the HIPASS survey (Koribalski et al. 2004) 
taken from a pointing at 16h47m20.0s, -57°27’06” (J2000). The y-axis is 
flux in Jy, and the x-axis heliocentric velocity in kms“^. The Hi velocity is 
842kms“^ taken at the midpoint of the Hi profile at the 20% level which 
has a velocity width of 222 km s“ ^. 
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Figure 9. The distribution of velocities of the star forming regions that are 
located along the ring (black) and further out or internal to the ring, includ¬ 
ing from component B and C (red). Circles are WiFeS data and triangles 
DBS. The horizontal axis shows the angle along the ring, starting from 
the direction of the assumed bullet, galaxy componentB.On the RHS the 
HIPASS spectrum has been plotted for ease of comparison. 


assumed bullet) including for knots labelled I and II in Fig|^ The 
average of these is 805kms“^, with little dispersion and is sepa¬ 
rated from that of the ring at this position angle of ~830km s“^. 


5.4 Masses 

5.4.1 Derived stellar mass 

The system’s stellar mass can be derived from the WISE drizzle 
images using the relation of Cluver et al. (2014), calibrated from 
the stellar masses in the GAMA survey (Taylor et al. 2011). The 
relation is: 

log ^ =-2.54 (lEl - IE2)-0.17 (1) 

Lwi 

where Lvvi is the in-band luminosity, i.e. the luminosity of the 
source as measured relative to the Sun in the W1 hand as defined 
in equation (1) in Cluver et al. (2014). 

Using Lwi = 8.27 X 10^ L©, we obtain a stellar mass of 
M* = 2.4 X 10® M© for the entire system. The colour correction 
used in this method aims at accounting for the contribution from 
a young stellar population. Without this correction derived masses 
would have been a few times higher. 

The main galaxy component stellar masses can also be esti¬ 
mated from their TTs-hand luminosities. The Jfs-band magnitudes 
were measured using 2MASS aperture photometry, subtracting the 
field stars within the aperture. The it's-hand image (Fig|^ clearly 
shows the main galaxy and a weak hut definite detection of galaxy 
components. The star forming ring is not discernible but evidence 
of a faint wash of near-IR flux internal to the ring is seen as well 
as a weak detection interior to the south-west edge of the ring as¬ 
sociated with component C previously identified as WKK 7457. As 


the TTs-band primarily shows the older stellar population, the main 
body of component A appears most prominent. 

Using the 2MASS Ks imaging and assuming a mass-to-light 
ratio of 0.80 (Bell et al. 2003), we find a stellar mass of the main 
components 1.29 x 10® Mq and 1.70 x 10® Mq for A and B, re¬ 
spectively. Their mass ratio is a factor of 7.5. 

The integrated TTs-magnitude of 8.95 gives a total stellar mass 
of 4.7 X 10® Mq, in broad agreement with the WISE derivation. 
The mass-to-light ratio at Ks varies by a factor of two over a wide 
range of stellar ages excluding the very youngest stars (Drory et 
al. 2004). The models of Henriques et al. (2011) show that the 
TTs-hand can be significantly enhanced by younger, more massive, 
AGB stars, and include a component from red supergiants. 


5.4.2 Gas and dust mass 


The Hi neutral hydrogen mass is reported by Koribalski et al. 
(2004) as Mhi = 2 x 10® Mq. This is similar to the overall stellar 
mass of the entire system which is therefore gas-rich. It is not clear 
where the bulk of the gas originates but the large Hi mass makes an 
original association with component A more likely. 

The mass in molecular gas is not known. CO data is required 
to determine this directly. Stark et al. (2013), provide evaluations of 
molecular to neutral Hydrogen (H 2/Hl) ratios from a hroad sam¬ 
ple of field galaxies spanning early-to-late types, evolutionary stage 
and stellar masses of 10^ ® to 10^® Mq. Their sample is not statis¬ 
tically representative because CO measurements for low-mass sys¬ 
tems such as ours are lacking. It does seem that for a range m g — r 
of ~0.3, molecular to neutral gas ratios range from 1 to ~10 per¬ 
cent so ~ 10 ® Mq in molecular hydrogen is possible. 

We use two different methods to estimate dust masses using 
available flux estimates. The AKARI 65, 90, 140 and 165p.m flux 
densities were measured using methods optimized for unresolved 
sources (Yamamura et al. 2009). Since the system is larger than 
the AKARI beam, the AKARI measurements may miss some of 
the extended emissioiQand could refer mainly to component A. 
This could explain why the AKARI flux estimates are lower by a 
factor of two compared to the equivalent IRAS values. We present 
AKARI fluxes for completeness but only use the WISE, IRAS, and 
Planck data in this analysis. 

The first method used fits the infrared SED with a blackbody 
function multiplied by an emissivity function that scales as . 
We then used the resulting temperature T in the expression: 


Mduat = 


K.B4T)’ 


( 2 ) 


where is a flux density measured in a single wave band, D is the 
distance to the source, is the wavelength-dependent opacity, and 
Bv{T) is the blackbody function. 

We perform fits to the data using j3 values of 1.5, as found em¬ 
pirically by Boselli et al. (2012) using Herschel data, and 2, which 
corresponds to the theoretical emissivity function given by Li & 
Draine (2001) that is frequently used in dust mass calculations. 
We also use the based on values from Li & Draine (2001) to 
calculate the dust mass corresponding to the modified blackbody 
fits with /3 = 2. Eor other /3 values it is common to rescale Kv 
by (A/Ao)~^, where Ao is some reference wavelength where the 


^ any contributions from the bright stars to the south are negligible, e.g. the 
SED for HD 150915 star drops by >3 orders of magnitude between 1 and 
10/rm. 
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emissivity is known. Though many authors select 250 /rm as a ref¬ 
erence wavelength, it is unclear if this is appropriate, as the choice 
is usually poorly justified. We only calculate dust masses for fits 
where j3 = 2, hut still report the temperatures for both /3 = 1.5 
and (3 = 2. The dust masses from our fits are derived from the am¬ 
plitudes of the best fitting modified blackbody functions. We use a 
function for where Ao=240 /rm and k„( 250 /im)=5.13 cm^ g“^. 

Fits were performed to the 100, 350, and 550p.m data that 
tends to originate from the coldest and most massive dust com¬ 
ponent heated by the diffuse interstellar radiation field. Emission 
at shorter wavelengths may come from a warmer thermal compo¬ 
nent (Bendo et al. 2006, 2012, 2015). The lOOtrm emission may 
also originate in part or entirely from dust heated by star formation. 
We need the data point to constrain the Wien side of the spectral 
energy SED, but we are also aware that the temperature may be 
overestimated by our analysis and the dust mass underestimated. 

We calculate the dust mass using the flux density at 350 /rm, 
the longest wavelength at which we have a > Scr detection. At 
shorter wavelengths, small uncertainties in the temperature will re¬ 
sult in relatively large uncertainties in the dust mass. We adopt 
K^{350fim) = 2.43 gcm“^ for the opacity. 

We obtain dust temperatures of 19 ± 2 K for fits with /3 = 2 
and 23 ± 3 K for /3 = 1.5. The slope of the 350-550p.m data is 
more consistent with a /3 = 2 emissivity function than the /3 = 1.5 
emissivity function found by Boselli et al. (2012), but given the 
large uncertainty in the 550tim data, it would be difficult to rule 
out the possibility that /3 is smaller than or larger than 2. The dust 
mass for the fit corresponding to /3 = 2 is (1.1 ± 0.3) x 10^ M©. 

In the second method, we fit the model templates of Draine 
et al. (2007) to the data. These templates are calculated for several 
different environments (the Milky Way, the EMC, and the SMC), 
different PAH mass fractions and different illuminating radiation 
fields. Draine et al. (2007) suggest fitting spectral energy distribu¬ 
tion data using the following equation which we also adopt: 

fv = A [(1 — 7)yiy((7MiN) + 7ii/([7MiN, (Tmax)] (3) 

In this equation, j,/([7min) is the template for dust heated by 
a single radiation field (7min and is equivalent to cirrus heated by 
a diffuse interstellar radiation field. The quantity ((7min , Dmax) 
is the template for dust heated by a range of illuminating radiation 
fields ranging from (7min to Umax, with the relation between the 
amount of dust dM heated by the a given radiation field in the range 
between U and U + dU given by the equation: 

dM oc dU (4) 

This yiy((7MiN, Dmax) effectively represents the fraction of 
dust heated by photodissociation regions (PDRs) in the fit and the 7 
term indicates the relative fraction of dust heated by the PDRs. The 
A term is a constant that scales the sum of the ji, functions (which 
are in units of Jy cm^ sr H-nucleon“^) to the flux density. It is 
related to the dust mass by: 

Mdust = (5) 

iVigas 

where Mdust/Mgas is the dust-to-gas ratio, mp is the mass of the 
proton, and D is the distance. 

As this system shares many characteristics of dwarf galaxies 
and Magellanic systems (including the metallicity - see below), we 
only fit EMC and SMC templates to the data. We use all WISE, 
IRAS, and Planck data between 12 and 550pm in the fit. 

We obtained best fits using templates based on the EMC2_05 
model from Draine et al. (2007). The resulting (7min and Umax 


were 1.5 and 10® times the local interstellar radiation field. The best 
fitting 7 was 0.024 ± 0.017. This is consistent with the majority 
or even all the dust mass being contained in a cirrus component. 
The resulting dust mass is (1.5 ± 0.3) x 10^ Mq, slightly higher 
than that calculated using the single modified blackbody fit. This is 
expected given that the Draine et al. (2007) models account for the 
mass in warmer dust components and for the relative contributions 
of various thermal components of dust to the SED. The masses 
are summarised in Table Individual stellar masses of the two 
systems are ~3 times less than that of the total system, so much of 
the stellar mass is located outside of the brighter regions. Deeper 
mid-IR images to find the detailed location of the remainder, either 
in a halo or in expelled streams, would be of interest. 

The two galaxy components A and B are small systems. For 
comparison, the stellar mass of the EMC is M*(LMC) = 2.7 x 
10® Mq. Its W1 flux is 2080 Jy (Melbourne & Boyer, 2013) which 
scaled to 10 Mpc gives 50mJy, about three times fainter than 
Kathryn’s Wheel. Component A is similar to the EMC. Compo¬ 
nent B is comparable to the SMC, which has a stellar mass of 

3.1 X 10® Mq. 

The system contains more gas and dust than the EMC and 
SMC. The EMC Hi mass is only 0.48 x 10® Mq (Staveley-Smith 
et al. 2003), and the SMC has an Hi mass of 5.6 x 10® Mq (Van 
der Marel, Kallivayalil & Besla 2009). While the stellar masses are 
broadly comparable, the dust mass is three times higher and the 
neutral gas four times higher than that of the EMC and SMC com¬ 
bined. We do not know with certainty which of the two components 
(A or B) is the origin of the gas and dust emission. The gas-to-dust 
ratio is ~133, which excludes the unknown molecular mass and is 
therefore an upper limit. For comparison, a mean gas-to-dust ra¬ 
tio of 58.2 was found for 78 bright galaxies in the Herschel Virgo 
Cluster Survey (Davies et al. 2012). Based on their study they sug¬ 
gest that a factor of ~ 2.9 can be used to convert atomic gas to total 
gas (i.e. Hi -I- H 2 ) in their survey. For comparison, the value of the 
gas-to- dust ratio for the Milky Way is ~143 (Draine et al. 2007), 
a similar value to what we derive here. Interestingly, the ratio of 
the extent of the star forming regions within B to that of the ring 
around A is similar to the ratio of the galaxy stellar masses. 

5.5 Star formation rates 

Star formation rates (SFR) can be derived from different tracers 
(e.g. Kennicutt & Evans 2012), mostly related to the UV photons 
produced by young massive stars. The Ha luminosity and the mid- 
IR flux tracing heated dust are independent tracers of the SFR. 

5.5.1 SFR from the Ha Flux 

The SFR from the integrated Ha luminosity of the system is calcu¬ 
lated following Kennicutt & Evans (2012), using the expression: 

log 1 = log - 41-27 ( 6 ) 

Mq yr-1 

Our first estimate of the integrated Ha flux (and luminosity) is 
derived from aperture photometry on all of the clearly identifiable 
knots in the CTIO continuum-subtracted Ha image. These included 
43 knots in the ring, 8 knots in the core of component B, and 8 knots 
in the region of the disk that defines component A. The resulting 
fluxes were converted to individual star formation rates for each 
knot. 

Aperture photometry was also done on the full areas covered 
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by each of the main system components A, B, C and the ring. This 
provided a global integrated SFR as Ha emission which is dis¬ 
tributed more widely around the entire region and even at low levels 
along the main disk of component A. The results are listed in Table 
For component B we measured the main body, but excluded the 
faint outer regions because of confusion with the ring to the south. 
Table 1^ also lists the specific SFR (SSFR: the SFR divided by the 
stellar mass). For the disk system component A and for B we used 
the stellar masses given in Table No mass was derived for the 
ring. For the full system we included all the stellar mass. 

The conversion factor of flux units count“'^ was derived from 
the individual Hll knots listed in Table using a 7 arcsec aper¬ 
ture for each. The Ha flux for these knots was measured from 
the calibrated WiFeS IFU data. Poor weather during the WiFeS 
observations resulted in a factor of two spread in measured con¬ 
version factors. Removing the lower values left knot systems I, 
II and IV. Averaging their conversion factors yields 1.2 ± 0.3 x 
erg cm~^ s“^ per count for the CTIO image. The photome¬ 
try counts were converted to flux units, and a total flux for all com¬ 
ponents of FfHa) = 3.0±0.6 X 10“^^ ergcm”^ s“^ was obtained. 

Owing to possible problems with the flux calibration, an inde¬ 
pendent measurement would be useful. Unfortunately, the system 
is just outside the footprint of the VST/OmegaCAM Photometric 
H-Alpha Survey: VPHAS-I-, (Drew et al. 2014), and while it is im¬ 
aged by the calibrated low-resolution the Southern H-Alpha Sky 
Survey Atlas (SHASSA) survey (Gaustad et al. 2001), the proxim¬ 
ity of HD 150915 precludes the measurement of an accurate inte¬ 
grated Ha flux following the recipe of Frew et al. (2013). How¬ 
ever, the SHS (Parker et al. 2005) offers an alternative approach to 
obtaining the integrated Ha flux of the system, in order to com¬ 
pare with our CTIO estimate. The SHS was absolutely calibrated 
by Frew et al. (2014a), who showed flux measurements to 20- 
50% accuracy were feasible for non-saturated diffuse emission re¬ 
gions. Using the methods discussed there, and correcting for re¬ 
duced transmission due to the systemic velocity, we determine an 
integrated Ha flux of 1.1 ± 0.5 x 10~^^ ergcm“^ s“^ for com¬ 
ponent B through an aperture of 90arcsec; a correction for [Nil] 
contamination is unnecessary based on the high excitation and low 
metallicity of the system’s Hit regions. This agrees with an equiv¬ 
alent CTIO measurement to within the uncertainties, so we aver¬ 
age the results to get a total flux for component B. Unfortunately 
the proximity of the bright star precludes an accurate flux mea¬ 
surement from the SHS for the star-forming ring. An approximate 
lower limit, based on a simple curve-of-growth analysis of several 
apertures is F’(Ha) = ~1 x 10“^^ ergcm“^ s”^, also in agreement 
with the CTIO measurement. The total Ha flux for the system, now 
using the weighted average of the independent CTIO and SHS cal¬ 
ibrations, is 2.9 ± 0.3 X 10“^^ ergcm“^ s“^, adopted hereafter. 

The systemic flux, now corrected for a reddening E{B—V) 
= 0.24 mag, and applying a Howarth (1983) reddening law, is 
Fo(Ha) = 5.1 X 10“'^^ erg cm“^ s“^. This leads to a total luminos¬ 
ity, L(Ha) = 6.1 X 10'*° ergs s“*, and a total ionizing flux, Q = 4.4 
xl0°° photons s“*, for a distance of lOMpc. The integrated SFR 
was then found to be 0.33 ± 0.09 Mq yr“* following equation|^ 
This SFR is about 60% higher than that of the LMC, adopting the 
Ha luminosity for the LMC from Kennicutt et al. (1995). Our es¬ 
timate may be an underestimate as we have not corrected for any 
foreground reddening, though the extinction derived from the mean 
Balmer decrement of the brightest Hll regions from our spectra is 
consistent with the Schafly & Finkbeiner (2011) extinction, indi¬ 
cating that the internal reddening is, in at least the most luminous 
Hll regions, minimal. There may also be photon leakage from the 


larger Hll regions in the system as noted by Relano et al. (2012), 
which might help to explain the difference between the SFRs de¬ 
rived from Ha and 22 qm 'WISE flux (see below). 

We note that the most luminous Hit region complex in 
ESO 179-13 is centred on Knot VII, with I/(Ha) ~ 5 x 10®° 
ergs s“*. This luminosity is comparable to supergiant HII regions 
like 30 Doradus in the LMC and NGC 604 in M33 (Kennicutt 
1984). Knot VII is within component C which contains old stars 
as it is detected in the 2MASS imagery as a separate object. It also 
lies interior to the fitted ring and has a different velocity distribu¬ 
tion to the components along the ring, exhibiting the largest veloc¬ 
ities across the system. We still consider it likely to be a separate 
star-forming dwarf galaxy but further observations are required to 
establish this clearly and to study the relationship between compo¬ 
nent A and C. 


5.5.2 SFR from the Mid-IR Flux 

The warm dust continuum is best traced by the 22 qm WISE flux. 
Cluver et al. (2014) have calibrated the relation of this flux against 
Ha-derived SFRs (corrected for extinction) from the GAMA Sur¬ 
vey (Gunawardhana et al. 2013) and find: 

M. 

log = 0.82 log vL{22 qm[L 0 ]) - 7.3. (7) 

This yields an integrated SFR of 0.47M©yr“* for ESO 179-13, 
somewhat higher than that derived from the de-reddened Ha flux. 
In contrast, the 12qm emission (Cluver et al. 2014), calibrated 
in the same way, corresponds to a lower SER of O.15M0yr“*. 
The lower mass measured in the 12 qm band may be because this 
band detects significant PAH emission. This is relatively weak in 
low metallicity environments (e.g. Engelbracht et al. 2005, 2008), 
so emission from bands that contain PAH emission and the SFRs 
based on these bands may be underestimated. All in all, the various 
estimates are in reasonable agreement, and we conclude that the 
integrated SFR is approximately 0.2-0.5 M© yr“*. 


5.5.3 Star formation distribution 

We find that almost all of the SFR in component A (the main disk 
or bar) is accounted for by the main identified knots. In contrast, 
only 28% of the SFR in component B appears concentrated in the 
main knots. About 44% of the SFR in the main ring is incorporated 
in the most intense and well defined knots. The remainder may be 
located in knots below the confusion limit of the available images. 

The SFR in the galactic disk that comprises component A is 
almost entirely due to the bright region towards the south west, lo¬ 
cated on the edge of the main disk and interior to the fitted oval 
of the main ring shown in Fig. There are a few further small, 
faint, reddened knots in the central regions but the remainder of 
the galaxy seems dormant. The star formation in component B is 
mainly concentrated in an oval region that is 30 x 20 arcsec in 
size that overlaps with the underlying galaxy detected in the optical 
broad band colours and 2MASS as galaxy WKK7463. However, 
there is star formation throughout a region 30-40 arcsec in radius 
centred on B with one small isolated blob even located 60 arcsec 
from the core to the north-west. There is clearly a continuum of 
star formation over much of the entire area where that associated 
with component B blends in to the north-east edge of the ring. 

The distribution of SFR over the detected knots is shown in 
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Table 5. Stellar, gas and dust masses, star formation rates and specific star formation rates estimated for the full system and, in so far as is known, for the 
individual components. Values for component A refer to the inner region (the disk), and for the core of componentB. The SSFRs for A and B assume total 
stellar' masses of 1.3 X 10® and 1.7 X 10® Mq respectively (refer to the text). 


Parameter 


Full system 

Component A 

Component B 

Ring 

Mt (Ks) 

[Mq] 

4.7 X 10® 

1.3 X 10® 

1.7 X 10® 


Mt (WISE) 

[M0] 

2.4 X 10® 




Mhi 

[Mq] 

2 X 10® 




Mduat 

[Mq] 

1.5 X 10^ 




SFR (Ha) 

[Mq yr-1] 

0.33 

0.01 

0.08 

0.24 

SFR (WISE 12 pm) 

[Mq yr-i] 

0.15 




SER (WISE 22 pm) 

[MQyr-i] 

0.47 




SSFR (Hq) 

[MQyr-iMT^l 

7.0 X 10““ 

7.7 X 10“^® 

4.7 X 10-1® 



Fig. The top panel shows the sum of the knots in compo¬ 
nents A, B and the ring, and the bottom panel shows the individ¬ 
ual components. The brightest knot in the ring (III.l) has a SFR of 
0.06 Mq yr“^. Galaxy componentB lacks the tail to low SFR seen 
in the knots of component A and in the ring: this is likely due to 
the confusion limit being much brighter for B. The two strongest 
regions of SFR are in the ring. 


5.5.4 Star forming versus non star forming galaxies 

Schiminovich et al. (2007) have shown that galaxies distribute over 
two clear sequences when the total stellar mass is plotted against 
specific star formation rate. The approximate location of the se¬ 
quences is shown in Fig. [m the dashed lines show the l-cr spread 
of the distribution. We plot the locations of component A compo¬ 
nent B and the total system. The SSFRs are taken from Table|^ The 
masses are taken as the full stellar mass divided over components A 
and B according to the ratio of their individual stellar masses. Note 
that component A includes the full mass but its SSFR only covers 
the apparent disk. 

ComponentB falls well above the star-forming sequence 
while component A fits close to the non-star forming galaxies as 
expected from a naive view of their star-forming activity evident in 
the imagery. Component A would be seen as quiescent if it were 
not for the bright region of star-formation at the south-west edge of 
the disk. The SSFR may be underestimated if extinction is signif¬ 
icant. Combining the entire system and using the full stellar mass 
and including the ring gives a data point at the lower end of the 



SFR [log M<,yr-‘] 


Figure 10. The distribution of star formation rates within the knots, for 
component A (along the disk), componentB (the assumed ’bullet’) and the 
ring. The top panel shows the combined distribution. 



Figure 11. The distribution of star-forming and non star-forming galaxies in 
the stellar mass versus specific star formation plot from Schiminovich et al. 
(2007) and adapted by Kennicutt & Evans (2012). Squares show locations 
of galaxies A, B and the full system. The ring is not included in the SFR for 
component A but is in the full system. 

star-forming sequence. The combined system is not vigorous, even 
though individual components in the system are. 

5.6 Abundances 

The observed optical emission-line spectra of the knots can be 
used to derive an approximate metallicity, using the oxygen abun¬ 
dance as a proxy. The weakness of the [N II] lines flanking Ha in 
each spectrum indicates low metallicity. Modelling of the spectra 
is made difficult by the lack of the [Out] 4363A line which is 
the main electron temperature indicator. However, metallicities can 
also be approximated from the [N II] and [S II] lines which are well 
detected. Fig.plots the line intensities relative to Ha of [Nil] 
and [S II], where [NII] refers to the 6584A line (the brightest of the 
doublet) while [S II] refers to the sum of both the 67I3A and 673lA 
lines. This plot, like the commonly used ‘BPT’ diagram (Baldwin, 
Phillips & Terlevich 1981), is used as a diagnostic to separate dif¬ 
ferent types of emission-line nebulae, but Viironen et al. (2007) 
and Frew & Parker (2010) have recently shown its sensitivity as a 
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metallicity and ionization parameter diagnostic. This diagram also 
circumvents the use of the [O III] line which is not always detected, 
and uses a small wavelength range in the red which improves the 
relative line-to-line calibration and extinction estimates. The [Nil] 
and [SII] lines are also produced at similar temperatures, alleviat¬ 
ing the sensitivity of the collisional forbidden lines to metallicity. 
The plotted diagnostic is a combination of the ‘N2’ and ‘S2’ indi¬ 
cators of Pilyugin, Vilchez & Thuan (2010). 

Fig.[^ shows the location of spectral line ratios from Hll re¬ 
gions measured across a variety of galaxies, in order of decreasing 
metallicity: M 31, the Milky Way, M 33, the LMC and lower metal¬ 
licity dwarf systems in the Local Group. The data are taken from 
Viironen et al. (2007, and references therein) and Frew et al. (in 
preparation). The boundary lines are adopted from Frew & Parker 
(2010) and Frew et al. (2014b), with the black-dashed trapezoid 
showing the domain of Hll regions. The plot shows a clear trend 
with metallicity with the different systems being well separated. 
Note that red curves enclose the planetary nebula domain (gener¬ 
ally of higher excitation) and the blue-dash enclosed area at the 
lower left is occupied by shock-excited supernova remnants (Sabin 
et al. 2013). The filled green diamonds show the location of the 
individual knots measured in ESO 179-13. Their location in the 
plot coincides well with Hll regions in the LMC, and therefore 
the galaxies can be considered to have a similar metallicity. In¬ 
deed, the relations of Pilyugin et al. (2010) indicate a metallicity 
of [O/H] ~ —0.4 for ESO 179-13, consistent with a value for the 
LMC. This makes it likely that the majority of the gas feeding the 
knots originated from galaxy A, which has a similar total mass and 
luminosity to the LMC. Components has a mass similar to the 
SMC and may therefore be expected to have a lower metallicity 
and gas content. 



0.0 ' . .. ' 

0.0 0.5 1.0 1.5 2.0 

log (Ha / [Sll]) 


Figure 12. Line ratios plotted for Hll regions in several galaxies: Milky 
Way (filled black squares), M31 (filled red squares), M33 (open red 
squares), LMC (filled blue squares), other LG dwarf galaxies (open blue 
squares). Refer to the text for further details. The ESO 179-13 emission 
knots are shown by the large green diamonds. 


5.7 Ages 

When did the bulTs-eye encounter occur between galaxies A and 
B? An age estimate can be derived from the projected angular sep¬ 
aration between the two components of 3.9 kpc and the observed 
heliocentric velocity difference between them. We assume they are 
effectively at the same distance and that the observed difference 
in the heliocentric velocities is due to their peculiar motions and 
that there is no relative acceleration. No transverse velocity com¬ 
ponent can be obtained so the estimated age is a lower limit. We 
use the average of the three literature value for component A from 
Strauss et al. (1992), Di Nella et al., (1997) and Woudt et al. (1999) 
of 761 kms“^ together with our newly determined value for com¬ 
ponent B from the average of several prominent emission knots in 
the system which gives a velocity of 806 ± 32kms“^ from fits to 
13 emission lines across the blue and red spectral regions. There is 
a nominal 45 km s“^ difference which produces an encounter age 
of 87 million years. If we assume the transverse velocity is of the 
same order as that along the line of sight then this introduces a root 
2 factor leading to a revised age of 123 million years. System age 
can also be estimated from the expansion assumed for the ring. Us¬ 
ing an assumed inclination of 45 degree, and a velocity amplitude 
(de-projected rotational/expansion velocity) of 85kms“^ then an 
age of 30 Myr is obtained. This is a lower limit as the velocities are 
mostly due to rotation, and the above assumes all is expansion. 


6 DISCUSSION 

6.1 Dwarf rings 

ESO 179-13 is among the smallest collisional ring galaxies known. 
The ring diameter of 6.2 kpc is the smallest of the Madore et al. 
(2009) catalogue, where diameters range up to 70 kpc. Compo¬ 
nents A and B are separated by only 4 kpc, while in the full cat¬ 
alogue, the colliding components are almost always more than 
lOkpc apart (Conn et al. 2011). The absolute K-band magnitude 
of ESO 179-13 is Mx = -21.1. Madore et al. (2009) list only 10 
objects that are less luminous. 

Because ESO 179-13 shows a full ring it is a unique case of 
a dwarf cartwheel-type galaxy. It is more than 1.5 mag (a factor 
of four) fainter than the next faintest example, Arp 147 (Mk = 
—22.8mag). It shows that the process causing full rings can also 
act in interacting dwarf galaxies. Such galaxies have smaller cross 
sections and shorter dynamical time scales, which may make such 
occurrences less frequent, but does not stop ring formation alto¬ 
gether. 

The difference in iT-band magnitude between components A 
and B is 2.2 mag. This is at the very top end of the distribution in 
Madore et al. (2009). As component A is already one of the small¬ 
est galaxies in the Madore et al. sample, component B is possibly 
the smallest of all known bullet galaxies causing (partial) rings. The 
mass ratio between A and B may be close to the limit for the for¬ 
mation for collisional rings formation (D’Onghia et al. 2008). 

6.2 Environment 

Collisional ring galaxies tend to be found in low-density environ¬ 
ments and avoid dense clusters (Madore et al. 2009). Kathryn’s 
Wheel as ESO 197-13 is also seen to be in a region devoid of ma¬ 
jor galaxies. The SHS images show no Ha-emitting companions 
within a degree of ESO 179-13. The broadband SuperCOSMOS R- 
band data reveal no other obvious galaxy of any kind within half a 
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degree of the centre of Kathryn’s wheel. There is a HIPASS galaxy 
with an almost identical velocity 1.5 degrees away: J1636-57 but 
is much fainter, with an integrated flux of 2.8 Jy.kms“^. Given the 
similarity in velocity and the angular separation of 260 kpc, a phys¬ 
ical association is possible. Galaxy J1620-60 also has the same ve¬ 
locity, is a few times fainter (37.4 Jy.km s“^) than ESO 179-13, but 
is over seven degrees, over 1 Mpc, away. The nearest large galaxy 
is NGC 6744, a grand-design barred spiral (Ryder, Walsh & Malin 
1999). It has an identical velocity and assumed distance (Kankare 
et al. 2014) and a projected separation of 19.1 deg, or 3.2 Mpc from 
ESO 179-13. It is unlikely to have interacted with this system. 


where both colliding galaxies formed rings though they are very 
different in nature. Table|^gives a summary table of the remaining 
key system parameters previously available or newly compiled. 

ESO 197-13 is an important system that can be used to study 
collisional rings as it is much closer than any equivalent, opening 
new possibilities for analysis. System dynamics can be studied at 
unprecedented spatial resolution as 0.1 arcsec (e.g. HST and/or AO 
resolution) is only 25 pc. This is sufficient to explore details of the 
system-wide star formation process and to allow detection of in¬ 
dividual stars. It would be a prime candidate to measure the star 
formation history in the regions recently traversed by the ring. 


6.3 Impact parameter 


The ring’s centre is displaced from that of of component A by about 
580pc. We take this as the impact parameter of the collision with 
the bullet components. The chances of a random (first) encounter 
between two isolated systems so close to the centre of A seems 
remote. Assuming three galaxies within a sphere of 250 kpc radius, 
and a relative velocity of 100kms“^, the chance of such a close 
encounter within a Hubble time is approximately lO”'^. Assuming 
2 X 10® galaxies within 10 Mpc, the chance of one such a close 
encounter happening within 10 Mpc within the past 10® yr (the life 
time of the ring) is 15%. This calculation suggests that ESO 197-13 
is likely the closest collisional ring galaxy to the Milky Way. 

The calculation ignores gravitational focussing. The gravity 
of the main galaxy bends the trajectory of the bullet towards it, in¬ 
creasing its capture cross section. This decreases the resultant im¬ 
pact parameter h. For v = lOOkms”^ and M — 6.6 x 10® Mq, 
the reduction is given by: 


b 

bo 



2GM\~^ 

~W) 


0.7 


( 8 ) 


This doubles the chances to 30%, and suggests a space density 
of 7 X 10“® Mpc“®. This estimate shows that the existence of one 
such system within 10 Mpc is quite plausible. 


7 SUMMARY 

During visual scans of the SHS (Parker et al. 2005), a spectacu¬ 
lar new collisional ring galaxy (or Cartwheel analogue) has been 
found, ESO 179-13, which we have named Kathryn’s Wheel. With 
an estimated distance of only 10.0 Mpc, it is in the local superclus¬ 
ter but has been previously missed due its location near the Galac¬ 
tic plane and its projected proximity to a bright 7.7-magnitude star 
which interferes with observation. 

The velocity range of all optical emission components from 
the main ring, component B and the published velocity for compo¬ 
nent A is about 200 km s“^. This is completely encapsulated within 
the FWHM of 222kms“^ reported for the Hi profile at the 20% 
level centred at a heliocentric velocity of 842kms“^ which we 
consider representative of the overall interacting system. 

At a distance of 10 Mpc, the system has an overall physical 
size of ~ 15 kpc for all major components and the diffuse surround¬ 
ing envelope, a total mass of M, = 6.6 x 10® Mq (stars -l- Hi) and 
a metallicity of [0/H]~ -0.4, that classifies it as a Magellanic-type 
system. It has a large reservoir of neutral gas and is the 60'’*' bright¬ 
est galaxy in the HIPASS Hi survey. By some margin it is also the 
nearest ring system known. This system offers an unprecedented 
opportunity to undertake detailed studies of this rare phenomenon, 
where two reasonably well matched galaxies have a close encounter 
centred close to the potential well of the primary. 


6.4 Lessons and prospects 

The discovery of Kathryn’s Wheel shows that the processes form¬ 
ing full rings can act in Magellanic type systems. Full rings have 
only been seen and studied in large, spiral galaxies, and so their 
formation has been assumed to require systems with halo masses 
above 5 x Mq (e.g. D’Onghia et al. 2008) though the cata¬ 
logue of Madore et al. (2009) does seem to extend to a few galaxies 
as small as ESO 197-13. As dwarf galaxies outnumber other Hub¬ 
ble types by a factor of ~ 1.4 (Driver 1999), dwarf collisional rings 
may also outnumber those systems involving spiral galaxies. This 
is shown by the space density derived above: it is 10 times higher 
than that of Few & Madore (1986), whose numbers are derived 
from more luminous galaxies. 

Both components appear to have been gas-rich prior to colli¬ 
sion. We assume the ring is formed in gas associated with compo¬ 
nent A, and so was gas-rich even though it shows little star forma¬ 
tion now (as is normal for residual galaxies within collisional ring 
systems). The collision may have diverted its gas out of the main 
disk. It is noteworthy that component B shows vigorous star forma¬ 
tion and may show its own (small) ring. It survived the encounter 
much better than A whose future evolution will be dictated by the 
depletion of its original gas reservoir even if the underlying stel¬ 
lar population seems largely unperturbed. Arp 147 is another case 
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Velocity Hi (helio) kms“^ 
Velocity (gal-centric) kms~^ 
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Distance (LG) Mpc 

842 ± 5 

736 ±6 

678 ± 11 
10.0 ±0.7 

Max. system extent (arcsec) 

318x198 

Max. system extent (kpc) 

15.4x9.6 

System metallicity [0/H] 

—0.4 

System age (Myr) 

30-120 

Space density (Mpc~^) 

7 x 10-5 







The nearest collisional ring galaxy 17 


Propulsion Laboratory/Califomia Institute of Technology, funded 
by the National Aeronautics and Space Administration, and data 
from the 2MASS, a joint project of the University of Massachusetts 
and the Infrared Processing and Analysis Center California Insti¬ 
tute of Technology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. The authors 
acknowledge help of Dr Anna Kovacevic at the CTIO telescope. 
We are grateful for stimulating discussion with Kate Zijlstra. We 
thank an anonymous referee who provided excellent feedback to 
help improve the quality of this paper. 


REFERENCES 

Amram P., Mended de Oliveira C., Boulesteix J., Balkowski C., 1998, 
A&A, 330, 881 

Appleton P.N., Stmck-Marcell C., 1996, Fund. Cosmic Phys., 16, 111 
Arp H.C., 1966, ApJ, 140, 1617 

Arp H.C., Madore B., 1987, A Catalogue of Southern Peculiar Galaxies and 
Associations. Cambridge, New York: Cambridge University Press 
Baldwin J.A., Phillips M.M., Terlevich R., 1981, PASP, 93, 5 
Barnes J.E., Hernquist L., 1992, ARA&A, 30, 705 

Beichman C.A., Neugebauer G., Habing H.J., Clegg P.E., Chester T.J., 
1988, The Infrared Astronomical Satellite (IRAS) Catalogs and At¬ 
lases. Volume 1: Explanatory supplement 
Bell E.F., McIntosh D.H., Katz N., Weinberg M.D., 2003, ApJS, 149, 289 
Bendo G.J. et al., 2006, ApJ, 652, 283 
Bendo G.J. et al., 2012, MNRAS, 419, 1833 
Bendo G.J. et al., 2015, MNRAS, 448, 135 
Bianchi L., 2014, Ap&SS, 354, 103B 
Block D.L. et al., 2006, Nature, 443, 832 
Boselli A. et al., 2012, A&A, 540, 54 

Boselli A., Fossati M., Gavazzi G., Ciesla L., Buat V., Boissier S., Hughes 
T, 2015, A&A, 579, 102 

Brosch N., Vaisanen R, Kniazev A.Y., Moisav A., 2015, MNRAS, in press, 
arXiv: 1506.00409 

Burbidge E.M., Burbidge G.R., Prendergast K.H., 1964, ApJ, 140, 1617 
ButaR., 1995, ApJS, 96, 39 

Childress M. J., Vogt F. P. A., Nielsen J., Sharp R. G., 2014, Ap&SS, 349, 
617 

Cluver M. et al., 2014, ApJ, 782, 90 
Condon J.J., 1992, ARA&A, 30, 575 

Conn B.C., Pasquali A., Pompei E., Lane R.R., Chene A-N., Smith R., 
Lewis G.E., 2011, ApJ, 741, 80 

Corwin H.G., de Vaucouleurs A., de Vaucouleurs G., 1985, Univ. Texas 
Monogr. Astron., 4, 1 

Davies J.I. et al., 2012, MNRAS, 419 3505 
Dennefeld M., Lausten S., Mateme J., 1979, A&A, 74, 123 
de Vaucouleurs G., de Vaucouleurs A., Corwin H.G., Buta R.J., Paturel G., 
Fouque R, 1991, Third Reference Catalogue of Bright Galaxies, Vol¬ 
ume III. New York: Springer 

di Nella H., Couch W.J., Parker Q.A., Paturel G., 1997, MNRAS, 287, 472 

D’Onghia E., Mapelli M., Moore B., 2008, MNRAS, 389, 1275 

Dopita M., et al., 2010, Ap&SS, 327, 245 

Doyle M.T., et al, 2005, MNRAS, 361, 34 

Draine B.T. et al., 2007, ApJ., 663, 866 

Drew J.E. et al., 2014, MNRAS, 440, 2036 

Driver S.R, 1999, ApJ, 526, L69 

Drory N., Bender R., Feulner G., Hopp, U., Maraston C., Snigula J., Hill, 
G.J., 2004, ApJ, 608, 742 
Engelbracht C., et al., 2005, ApJ, 628, L29 
Engelbracht C., et al., 2008, ApJ, 678, 804 
Lew J.M.A., Madore B.E, 1986, MNRAS, 222, 673 
Eosbury R.A.E., Hawarden T.G, 1977, MNRAS, 178, 473 
Ereeman K.C., de Vaucouleurs G., 1974, ApJ, 194, 569 
Erew D.J., Parker Q.A., 2010, RASA, 27, 129 
Erew D.J., Bojicic I.S., Parker Q.A., 2013, MNRAS, 431, 2 


Erew D.J., Bojicic I.S., Parker Q.A., Pierce M.J., Gunawardhana M.L.R, 
Reid W.A., 2014a, MNRAS, 440, 1080 
Erew D.J. et al., 2014b, MNRAS, 440, 1345 

Gaustad J.E., McCullough PR., Rosing W., VanBurenD., 2001, PASP, 113, 
1326 

Gerber R.A., Lamb S.A., Balsara D.S., 1992, ApJL, 399, L51 
Gil de Paz A., Madore B.E, 2003, ApJS, 147, 29 
Graham J., 1974, Observatory, 94, 290 
Gunawardhana M.L.R et al., 2013, MNRAS, 433, 2764 
Hambly N.C. et al., 2001a, MNRAS, 326, 1279 

Henriques B., Maraston C., Monaco R, Eontanot E, Menci N., De Lucia G., 
Tonini C., 2011, MNRAS, 415, 357 
Higdon J.L., Higdon S.J.U., Rand R.J., 2011, ApJ, 739, 97 
Hinshaw G. et al., 2013, ApJS, 208, 19 
Howarth I., 1983, MNRAS, 203, 301 
Ishihara D. et al., 2010, A&A, 514, 1 
James RA. et al., 2004, A&A, 414, 23 

Jarrett T.H., Chester T., Cutri R., Schneider S.E., Huchra J.R, 2003, AJ, 
125, 525 

Jarrett T.H. et al., 2012, AJ, 144, 68 
Kankare E. et al., 2014, A&A, 572, A75 
Kennicutt R.C., 1984, ApJ, 287, 116 

Kennicutt R.C., Bresolin E, Bomans D.J., Bothun G.D., Thompson I.B., 
1995, AJ, 109, 594 

Kennicutt R.C., Lee J.C., Eunes S.J., Jose G., Sakai S., Akiyama S., 2008, 
ApJS, 178, 247 

Kennicutt R.C., Evans N.J., 2012, ARA&A, 50, 531 
Koribalski B.S. et al., 2004, AJ, 1, 16 

Lauberts A., 1982, ESOAJppsala survey of the ESO(B) atlas. Garching: 
European Southern Observatory 

Laustsen S., Madsen C., West R.M., 1987, Exploring the Southern Sky. 

Garching: European Southern Observatory 
Lee, J.C. et al. 2011, ApJS, 192, 33 
Li A., Draine B.T., 2001, ApJ, 554, 778 
Lindsay E.M., Shapley H., 1960, Observatory, 80, 223 
Longmore A.J. et al., 1982, MNRAS, 200, 325 
Madore B.E, Nelson E., Petrillo K., 2009, ApJS, 181, 572 
Marston A.R, Appleton P.N., 1995, AJ, 109, 1002 
Masci E, 2013, Astrophysics Source Code Library, 2010 
Masci F.J., Fowler J.W., 2009, Astronomical Data Analysis Software and 
Systems XVIII, 411, 67 

McDonald I., Zijlstra A.A., Boyer M.L., 2012, MNRAS, 427, 343 
Melbourne J., Boyer M.L., 2013, ApJ, 764, 30 
Meyer M.J. et al., 2004, MNRAS, 350, 1195 

Murphy T., Mauch T, Green A., Hunstead R.W, Piestrzynska B., Kels A.R, 
Sztajer R, 2007, MNRAS, 382, 382 
Parker Q.A., Bland-Hawthorn J., 1998, RASA, 15, 33 
Parker Q.A. et al., 2005, MNRAS, 362, 689 
Parker Q.A. et al., 2006, MNRAS, 372, 79 
Pilyugin L.S., Vilchez J.M., Thuan, T.X., 2010, ApJ 720, 1738 
Planck Collaboration: Ade P.A.R. et al., 2014, A&A, 571, 28 
Relano M., Kennicutt R.C., Eldridge J.J., Lee J.C., Verley S., 2012, MN¬ 
RAS, 423,2933 

Rodgers A.W., Harding R, Bloxham G., Bessell M.S., 1998, PASP, 100, 
626 

Romano R., Mayya Y.D., Vorobyov E.I., 2008, AJ, 136, 1259 

Ryder S.D., Walsh W., Malin D., 1999, RASA, 16, 84 

Sabin L. et al., 2013, MNRAS, 431, 279 

Schiminovich D., Wyder T.K., Martin D., 2007, ApJS, 173, 315 

Schlafly E.F., Finkbeiner D.R, 2011, ApJ, 737, 103 

Schroder A., Kraan-Korteweg R.C., Henning P.A., 2009, A&A, 505, 1049 

Skrutskie M.F. et al., 2006, AJ, 131, 1163 

Sersic J.L., 1974, Astrophys. Sp. Sci., 28, 365 

Smith R.T., 1941, PASP, 53, 187 

Smith R., Lane R.R., Conn B.C., Fellhauer M., 2012, MNRAS, 423, 543 
Stark D.V, Kannappan S.J., Wei L.H., Baker A.J., Leroy A.K., Eckert K.D., 
Vogel S.N., 2013, ApJ, 769, 82 


18 Parker et al. 


Staveley-Smith L., Kim S., Calabretta M.R., Haynes R.F., Kesteven M.J., 
2003, MNRAS, 339, 87 

Strauss M.A., Huchra J.R, David M., et al., 1992, ApJS, 92, 29 
Taylor E.N. et al., 2011, MNRAS, 418,1587 
Theys J.C., Spiegel E.A., 1976, ApJ, 208, 650 
Theys J.C., Spiegel E.A., 1977, ApJ, 212, 616 

Tully R.B., 1988, Nearby Galaxies Catalog. Cambridge and New York: 
Cambridge University Press 

van der Marel R., Kallivayalil N., Besla G., 2009, in Proc. lAU Symp. 256, 
The Magellanic System: Stars, Gas, and Galaxies, J.T. van Loon & J.M. 
Oliveira, eds., p 81 

Viironen K., Delgado-Inglada G., Mampaso A., Magrini L., Corradi 
R.L.M., 2007, MNRAS, 381, 1719 

Vogt F.P.A., Dopita M.A., Borthakur S., Verdes-Montenegro L., Heck¬ 
man T.M., Yun M.S., Chambers K.C., 2015, MNRAS in press 
larXiv:1504.03337l 

Wright A.E., Griffith M.R., Burke B.F., Ekers R.D., 1994, ApJ. Suppl. Ser. 
91, 111 

Woudt P.A., BCi'aan-Korteweg R.C., Fairall A.P, 1999, A&A, 352, 39 
Woudt PA, Kraan-Korteweg R.C., 2001, A&A, 380, 441 
Wright E.l. et al., 2010, AJ, 140, 1868 
Yamamura I. et al., 2009, ASP Conf. Ser., 418, 3 

Zwicky F, 1941, in Theodore von Karman Anniversary Volume: Contribu¬ 
tions to Applied Mechanics and Related Subjects. Pasadena: California 
Institute of Technology, p. 137 


